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Abstract. We study a particular type of local quench in a generic quantum critical one-dimensional 
system, using conformal field theory (GET) techniques, and providing numerical checks of the results in 
free fermion systems. The system is initially cut into two subsystems A and B which are glued together 
at time t = 0. We study the entanglement entropy (EE) between the two parts A and B, using previous 
results by Galabrese and Gardy ^[2], and further extending them. We also study in detail the (logarithmic) 
^ Loschmidt echo (LLE). For finite-size systems both quantities turn out to be (almost) periodic in the scaling 
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limit, and exhibit striking light-cone effects. While these two quantities behave similarly immediately after 
^ the quench — namely as c/31og4 for the EE and c/41og4 for the LLE — , we observe some discrepancy once 

ethe excitations emitted by the quench bounce on the boundary and evolve within the same subsystem A (or 
B). The decay of the EE is then non-universal, as noticed by Eisler and Peschel [3- In contrast, we find 
^ I ^ that the evolution of the LLE is less sensitive than the EE to non-universal details of the model, and is still 

accurately described by our CFT prediction. To further probe these light-cone effects, we also introduce a 
variant of the Loschmidt echo specifically constructed to detect the excitations emitted just after the quench. 
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1. Introduction 

The characterization of entanglement in quantum many-body systems has become increasingly important 
over the last few years. It has been at the heart of a fruitful interdisciplinary work, with strong motivations 
coming both from quantum information |3] and from "conventional" condensed matter theory. These include 
topics as diverse as the investigation of computational complexity of algorithms such as DMRG [S] [S] , or the 
development of new numerical tools such as PEPS jTJ or MERA [8 — which aim to give access to far larger 
system sizes than exact diagonalisation techniques — , or the hunt for theoretical tools able to grasp some of 
the features of topologically ordered phases of matter jU] [TO] [H] ■ 

Among all the competing quantities that have been proposed as (theoretical) measures of quantum 
entanglement, the von Neumann entanglement entropy (EE) is the most widely accepted [U [TH [T3]. It 
is defined through a bipartition of the Hilbert space T-L = Ha ®T-Lb- Most of the time, the bipartition 
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corresponds to a spatial cut of the system, the subsystem A being thought of as the physical system of 
interest, and B as the "environment". Usually, one starts with a pure state \ ^ H, then defines its reduced 
density matrix pA = Tr-^^ |^^) {^\, and its EE as the von Neumann entropy of pA 

S = -Tr^^p^logpA- (1) 

Intuitively, the more entangled the system is, the less separable it is, and therefore the larger the EE is. While 
the EE has been computed numerically or analytically in a large variety of systems, its experimental measure 
would generically require the knowledge of a prohibitively large number of correlation functions, which can 
be traced back to the intrinsic non-locality of this quantity. However, some methods have been proposed for 
direct measures in free fermion systems [T^ [T51 [TB] (see also the discussion in [TT] about interacting systems), 
or for an original experimental implementation of the replica trick which could possibly lead to a measure JS] 
of the Renyi entropies Sn = logTrp^. 

In [19) we introduced a quantity which has a behaviour that is similar to the EE, the (logarithmic) 
bipartite fidelity (LBF). The LBF is in principle simpler to measure than the EE, because it is defined as 
an overlap of two states (see below). The common feature of the "experimental" methods that have been 
proposed to measure both the EE and the LBF is that they rely on probing some of the non-equilibrium 
properties of the system. More specifically, the relation between entanglement and time-evolution of quantum 
systems after a sudden change of parameters in the Hamiltonian, i.e. after a quantum quench, has been 
emphasized. In that sense, the methods proposed for a measure of the Renyi entropy T8J and for the LBF 
[T9j are the same, and rely on the probability of observing the system in the ground state of the Hamiltonian 
H' after a sudden local quench H ^ H' . An experimental device would be considerably more involved in 
the case of the Renyi entropies [18] than for the LBF though. The theoretical study of quantum quenches 
has been very active [20l[2ll[22l[23l|24l[25l[26l[23|28] (the list is far from exhaustive), in part triggered by 
the recent experimental achievements in ultracold atoms [351 130] : which now allow to build atomic systems 
without dissipation, and with large ranges of tunable parameters. Another fascinating question is the one of 
thermalisation in quantum systems after quenches [3TJ [351 [Ml ISi] . 

One of the most interesting aspect of some of the results on entanglement and quantum quenches is their 
universality |36| , I37| , [38]. Various scaling laws, critical exponents and universal functions have been discovered. 
In the case of one-dimensional systems, in particular, it is well-known that, at quantum criticality, the low- 
energy excitations are generically described by a conformal field theory (GET). The GET formalism has been 
successfully applied to derive a variety of results in that context (see for example [J (THl HZ] and references 
therein) . 

In this paper, we fill some of the gaps that still appear in this already large literature. Indeed, to our 
knowledge, the question of finite-size systems and quantum recursion at criticality has not been addressed 
in general (see, however, [SOj for partial results). Neither has been studied the LBF after a local quench in 
such systems, apart from [19] . Therefore we focus on local quantum quenches in finite-size systems, which 
we study by GET methods, extending the results of the pioneering papers by Galabrese and Cardy P] and 
Eisler and Peschel [S] I41j and providing numerical checks of our results in free fermion systems. We study 
both the EE and the LBF, emphasizing the differences between the behaviours of these two quantities. We 
find that there are striking light-cone effects, both for the EE and for the LBF. Since for a system evolving 
in time the LBF is nothing but (minus the logarithm of) the Loschmidt echo 

c{t) = mt)\m)\\ (2) 

as noted in [19j . in the rest of this paper we will prefer the latter terminology, as it is more often used in 
the context of quantum quenches. We will therefore call logarithmic Loschmidt echo (LLE) the quantity 
F{t) = — log£(t), keeping in mind that it is the same as the LBF defined in ,19 evolving with time t. 

The paper is organized as follows. In Sec.|2]we first describe some basic consequences of a local quantum 
quench in a one-dimensional critical system. In particular, we explain why the time evolution of the wave 
function after the quench should be periodic. In Sec. |3] we derive several results for the time-evolution of 
the entanglement entropy as well as the Loschmidt echo. This section relies mainly on GET methods, and 
may be skipped at first reading. Finally, Sec. |4] focuses on the results, providing physical interpretations and 
numerical checks of these in free fermion systems. 
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2. Finite-size critical systems in one dimension 

In this paper we consider 1+1-dimensional systems which are described by a CFT at sufficiently large 
distances and time scales (see for example [42l [43] for an introduction to CFT) . We are interested mainly in 
finite-size systems. Typically, such a system can be a gapless spin chain of finite size L, for sufficiently large 
L. For a system with open boundary conditions, which are assumed to correspond to a single conformal 
boundary condition throughout the paper, the lower-energy part of the spectrum of the gapless Hamiltonian 
Hl is given by [HI SS] 

Ea = Li^bulk + -^boundary + {ha — — ^ + O (f/i^) (3) 

L i?buik is the (bulk) extensive part of the energy, the boundary energy. i?buik and Eboundary 

are both non-universal, vp is the Fermi velocity, c is the (universal) central charge of the CFT and the ha 
are the conformal dimensions of the operators (pa appearing in the low-energy part of the spectrum. 

There are two kinds of such operators, the primary fields and their descendants. The descendants can 
always be constructed out of a single primary field, and correspond to some combinations of its derivatives. 
The primary fields cannot in general be expressed as descendants of other fields. The primary fields that 
appear in the spectrum ([s]) depend both on the CFT and on the boundary condition we are dealing with. 
Typically, in unitary CFTs that are of interest for condensed matter applications (for example as scaling 
limits of spin chains), there is a discrete (often finite) set of primary fields. If a primary field is present in 
the spectrum (|3|, then all its descendants must appear as well. The conformal dimensions of a descendant 
4)a,k and of the primary field (pa are always separated by an integer. This property is crucial when we deal 
with time evolution. Indeed, if the system is initially in the state \ip{0)), then \ilj{t)) = e~*'^* \ip{0)), and 
in particular after time t = 2L/vf all the descendants of a primary field (pa pick up the same phase, up to 
additive 0(1/L^) x L = 0{1/L) corrections. 
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Figure 1. The "quasi-particle" picture gives a useful though heuristic way of thinking about the time- 
evolution of finite-size critical systems in Id. All the quasi-particles move at speed vp. After time 2L/vp, 
the particles come back at the same position. Therefore, up to possible phase factors due to the reflection 
against the boundaries, the system is periodic. 

This can be interpreted intuitively with the following semi-classical picture (see also Fig. [T]). Low-energy 
excitations can be seen as some "quasi-particles" moving at fixed speed vp- There can be different species 
a of quasi-particles in the system, which correspond to the different primary fields pa in the CFT. The 
behaviour of these different species of quasi-particles differs at the boundary: they pick up a phase factor 
g-i7r/ta^ This "quasi-particle" picture is heuristic and is useful to predict light-cone behaviours, as explained 
below, however it fails to give quantitative predictions in general. 
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It is important at this stage to emphasize that the commensurabiUty of the spectrum (more precisely 
of each subsector of the spectrum corresponding to each primary field (j)a and its descendants) can only 
occur in the scaling limit. It is not true in a finite system, because of the order 0{1/L^) corrections in (|3|. 
The crucial requirement is that the dispersion relation must be linear — vp\k\ (if k is the momentum 
in a translation- invariant system). In our numerical checks we will focus on finite-size lattice systems, and 
although the dispersion relation can be linearised close to the Fermi energy, there are subleading terms which 
spoil the commensurability of the spectrum (see [35] for a detailed discussion of this point). However, as we 
shall see, the periodic behaviour will be a very good approximation in the case of the local quenches we are 
interested in: see already Fig. [2] for a first example. There are two reasons why this is the case. 

• We look at the behaviour of the system at sufficiently large distance and time scales (typically large 
compared to the lattice spacing a), such that the continuum description is relevant. However, we restrict 
our analysis to times t ^ L'^ /(avp), such that the order 0{1/L^) terms in ^ are still small. For larger 
times, these terms spoil the periodicity. Typically, for systems with a number of sites of order L/a 100, 
we expect to see about 10 oscillations that are well described by the CFT calculations, and after that 
the non-universal corrections become too large. 

• We are interested in local quenches only. The Hamiltonian i?^ — > H'j^ is changed only locally. This 
is radically different from a global quench, where parameters of the Hamiltonian would be changed 
suddenly everywhere in the system. After a local quench, there are much less excitations in the system 
than after a global one. In particular, for a local quench, we expect only the lower energy part of the 
spectrum to be significantly excited. For these excitations the formula ([3| applies. For global quenches 
however, we would expect also high-energy excitations with non-universal scaling behaviour. Lattice 
effects would be presumably much more important in that case, and would spoil the commensurability 
of (the different sectors of) the spectrum almost immediately, leading to chaotic or turbulent behaviour, 
as observed in [Ml. 
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Figure 2. Entanglement entropy and the logarithmic Loschmidt Echo against rescaled time {vp/L)t after 
a local quench in the middle of a system of size L = 128. A couple of oscillations of the finite-size system 
can be observed after the local quench. They are very well described by the CFT (slightly anticipating 
on Eq. |33| and Eq. |37| . The rapid oscillations appearing in the plots are analogous to the one observed at 
equilibrium (see for example 1461 1471 l48l V Their amplitude decreases for bigger and bigger system sizes, and 
we will not discuss them in this paper. After longer times, the non-universal 0(1/-L^) terms in (jsjl spoil the 
CFT oscillations, and the periodic behaviour is lost. 



3. Local quenches from conformal field theory 

We now specify the type of local quenches we are interested in, and we show how these can be studied using 
CFT techniques, following previous works [21 HTl [T^ . 
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Figure 3. The system is initially prepared in the ground state of Hj^ + Hg, then at time t > it evolves 
with the total Hamiltonian H = Hj\ + Hg + H\„ . 



3.1. Cut and glue: the type of local quench 

In this paper we study one type of local quench only. Let us consider a one-dimensional system of length L, 
with open boundary conditions. We cut the system into two parts A and S, of length La and Lb = L — La 
respectively. The system is initially prepared in the ground state of the Hamiltonian Ha + Hb , which is of 
course the tensor product jA) ® of the ground state \A) of Ha and the ground state \B) of Hb- For later 
convenience we note \A ® B) = \A) ® \B) . At time t = 0, the interaction H\g between A and B is switched 
on instantaneously. The system then evolves with the total Hamiltonian H = Ha + Hb + H\q (see Fig.|3]). 

Throughout this paper we consider situations when the boundary condition (b.c) is always the same. 
Namely, if we are dealing with a spin chain of size L with free b.c, then the boundary conditions at the 
cut between A and B are also free. Moreover, we assume that this is a conformal b.c. This is always 

true at least in the scaling limit, since every b.c should renormalize towards a conformal b.c under the 
renormalization group (RG) flow. In this paper we address neither the question of different b.c, which would 
involve boundary condition changing operators, nor the one of boundary RG flows. 

There is an important consequence to these assumptions. In general, the initial state \A (g) B) should be 
decomposed as 



where the index a runs over the different primary fields in the spectrum of the GET, and k runs over the 
descendants of each primary field. This, in principle, gives the time-evolution after the quench, since by 
definition the \4>a,k) are the eigenstates of the Hamiltonian H, with eigenvalues ([3|. The calculation of 
the coefficients Ca^k is difficult though. For the particular type of quench we are looking at, however, only 
the identity and its descendants have non-zero overlaps Ca.k = {4'a,k\A<Si B). This huge simplification occurs 
because we assume that the conformal b.c. is the same everywhere. Since the initial state \AiSi B) = \ A)(Ei\B) 
is built out of the ground states of Ha and Hb, which are both Hamiltonians with the same b.c. on the 
left and right boundary, \A) and \B) both correspond to the identity operator. Then the fusion of these two 
states has to be the identity as well. This, composed with some conformal transformation of the space-time 
geometry (as will become clear below) cannot generate any operator but the identity and its descendants. 

Since only the sector of the identity appears in the spectrum, the evolution of the system should be truly 
periodic (up to the order 0(1/L^) corrections discussed in the previous section). There is no additional phase 
corresponding to different primary operators which could presumably affect the periodicity anymore. In the 
"quasi-particle" picture this means that there is only one species of quasi-particles which can be emitted by 
the local quenches we are considering. 

Let us stress again the fact that this conclusion would not hold if we were considering different b.c. in 
the system. Ghanging the b.c, for example on both sides of the quench, would generate other species of 
quasi-particles. 
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3.2. Space-time geometry and conformal mappings 



As a warm-up exercise, let us consider the simpler case of two semi-infinite systems A and B that are glued 
together at time t = 0. This of course corresponds to the limiting case a -€1 vp t La, Lb in the above 
setting (Fig. [s]) . We review the CFT calculation of the Loschmidt echo [T5] and of the entanglement entropy 
[2] in that case. 




Figure 4. The geometry used to define the partition function Za b- It can be obtained by a conformal 
mapping from the upper half-plane z i— >■ w(z) = (r/4)(2 -|- z~^). Notice that the imaginary time flows 
in the left-right direction. The distance between the two slits is r. The inflnitesimal transformation 
w I— >■ w{l + 5t/t) is applied at the right of T in the w region. F' is the image of F through the inverse 
conformal mapping: F' = ui~^(F). 



3.2. L Loschmidt echo The CFT calculation is performed in imaginary time r. For the Loschmidt echo we 
need to compute 

C{T)^\{A(g,B\e-^'' \A(g)B)\\ (5) 
where H = Haub is the Hamiltonian of the total system. This can be done using the relation 

\A(^B) = lim e-^(^-*+-f^«)e^-^° |s) / {s\A <E) B) , (6) 

A— f oo 

which is valid for any state \s) such that {s\A ® B) 7^ 0. is the ground state energy of Ha + Hb- The 
relation ([S]) becomes (for r > 0) 

2 



C{t) = est X 



lim Za,b{t) 



(7) 



where Za,b{t) = (s| e-^('f^-4+HB)g-rHg-A(ffA-(-HB) This is the partition function of a statistical 

mechanical system defined in the geometry shown in Fig. |4j The normalization constant in Q will be 
fixed later. 

There is an extensive part in the free energy J^a,b{t) = — log Za,b{t) coming from the region of size r 
between the two corners in Fig. |4] This part will contribute only to a phase factor when we come back to 
real time (r — )■ it), which disappears in the norm ([s]). Also, there is an extensive free energy which grows 
with the size of the system and diverges as A — t- 00. We can absorb this contribution in the normalization 
constant however, which will be fixed later when we come back in real time. Thus, the interesting part in 
the calculation of the time-evolution of the Loschmidt echo is precisely the universal part coming from the 
corners in Fig. |4] As already emphasized in [T^ , the simplest way to get the result at this point is to apply 
the Cardy-Peschel formula [SO]. However, we wish to compute the Loschmidt echo in more complicated 
geometries later on, and for pedagogical purposes we will follow a slightly different route. 

The dependence of J^{t) ~ — log£(r) on the distance r in Fig. |4] can be evaluated as follows 50 . 
Consider the transformation w ^ w + 5t ioT: w aX the right of the path T in Fig. |4j and w ^ w dX the left 
of r. Notice that this mapping changes the distance t into r + St. The variation of T can be expressed in 
terms of the stress-tensor T(w) (and its anti-holomorphic part T{w)) as 

= g / {T{w)) dw + c.c. (8) 



r 
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and the one-point function of the stress-tensor can be obtained from a conformal mapping from the half-plane 
(see Fig. l4| 



wiz) ^ -{z + z 1). 



(9) 



The one-point correlation function of {T{z)) is zero in the half-plane. {T{w)) can be calculated using the 
transformation law of the stress-tensor [421 143] 

c 



T{w) — (dw/dz) T{z) — S{w,z) 



12 

where S{w, z) is the Schwarzian derivative of the conformal transformation w{z) 
S{w, z) 



w"'{z) 
w'{z) 



w"{z) 
w'{z) 



We therefore get 



c6t 



dz S(w, z) 



dw\ 
dz I 



(10) 



(11) 



(12) 



where the path V is the image of F under the inverse map w ^ z. For the mapping in Eq. [9] the integrand 



f[z) = S{w,z) 



dw 
dz 



24 



(13) 



T (Z- l)3(z+l)3- 

Notice that f{z) has singularities at z = ±1, which are the images of the two tips of the slits in Fig. [4] through 
the inverse mapping w z. In this case we have F = and F' — The integral can be evaluated as 

follows. First we consider both F' = i]R+ and its mirror image: iM+ iM., which multiplies the integral by 
2. Then we close the contour by a big half-circle to the right. We call this contour C, then 



f{z) dz^ - 



f{z)dz = - 



f{z)dz. 



(14) 



Finally the integral over C can be computed using the residue theorem. Only the pole at z = 1 will contribute, 
and we obtain 



5T 

1^ 
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Res[/(z);z = l] 



(15) 



This equation is particularly useful because the computation of the Loschmidt echo for all the other 
geometries will lead to a similar result. Namely, the calculation boils down to a residue at the image of 
one of the tips through the inverse mapping. Here we get 6J^/6t = which gives finally 



^(r) = -log£(r) = -log|r| 



est. 



(16) 



Now we can perform the Wick rotation r — > ivpt + e and come back in real time. An UV cutoff e is needed 
here in order for the Loschmidt echo to be defined at time i = (note that £(t = 0) = 1). It is of the order 
of a lattice spacing in the case of spin chains. The logarithmic Loschmidt echo (LLE) is then 

F{t) = - log L{t) = ^ log Vl + {vpty/e'^ ^ ^ log i + est (17) 

when which is the result of [TO] . 



3.2.2. Entanglement entropy We recall quickly the calculation of the entanglement entropy, following 
Calabrese and Cardy[2l. The Renyi entropy 

^n = 7^1og(Trp"), (18) 
1 — n 

gives the EE in the limit n ^ \, assuming that Sn is analytic in n 

5 = -Trplogp= YixaSn (19) 

When n > 2 is an integer, Trp" admits a simple interpretation in terms of partition functions on a n— sheeted 
Riemann surface, and can be computed in CFT[T1 While the validity of the replica trick here is by no 
means obvious, the method has enjoyed considerable success so far and has proved to be very powerful in 
the context of ID quantum critical systems. 
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As explained in detail in [D H] , Tr transforms as a one-point function of a primary field $„ under a 
conformal transformation. Its conformal dimension is 



12 



1 



(20) 




Figure 5. The geometry used to compute the entanglement entropy. <3?n is tlie twist operator introduced 
by Calabrese and Cardy in [l]. 

The conformal mapping needed is the same as for the Loschmidt echo, with r changed into 2e (see Fig. [5]) 

w{z)^'-{z + z-'). (21) 

The Renyi entropy is given in this CFT calculation by 
1 



Sn 



1 — ri 



log ($„(w)) 



The one point function can be obtained by mapping back to the z— half plane (Fig. [5]) 

($„(w)) cx 



5m(z) dw 
and the entropy is expressed as 



1 dz 



(22) 



(23) 



(24) 



3m(z) dw 

Since the initial cut coincides with the bipartition, the operator lies between the two slits at some point 



w = 9 on the real axis, with — e < 6* < e (see Fig. [5|). Inverting (21 ), we get 
1 dz 



9m(z) dw 

Performing the Wick rotation 6 ivpt we recover the celebrated result [2] for the EE {n — 1) 



S{t) = -logt + cst 



(25) 



(26) 



in the limit t 3> e. 



3. 3. A simple finite size example 

3.3.1. Loschmidt echo We consider now the case of a cut in the middle of the chain La — Lb — L/2. The 
conformal transformation which maps the upper-half plane onto the "double-pants" geometry in Fig. [6] is 



w{z) 



L 
2^ 



log 



1 



z + l 



log 



(27) 



^z + a ^ 

This Schwarz-Christoffel mapping has singularities at z = ±l,z = ±a, and extrema at z = i\/a. The 
distance between the two slits is 

4L 

T— \w{\/a) w{—\/a)\ — — arctanh (l/-v/a) (28) 
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—a —y/a —1 1 y/a a 



Figure 6. Pants geometry (left) and upper half-plane H (right). The conformal transformation w{z) is 
given by Eq. [27] 



Inverting this relation, we get 
air) 

Denoting as before f{z) — S{w, z){dw/dz)^^ , the residue we have to calculate is then 
Res[/(z); z = ^a] 



Stt + 6a + 1 



8L y/a{a + 1) ' 

and the variation of the (logarithmic) Loschmidt echo with the imaginary time r is 

S log C TTc + 6a + 1 ttc 1 



St 6t 
Integrating over r, we get 



16L ^(a + 1) 4Ltanh(7rr/L) 



J-(r) 



log£(T) — - log |sinh(7rT/L)| + est. 



(29) 



(30) 



(31) 



(32) 



Of course, the constant in (|32| can still have a non-trivial dependence in L. This is indeed the case as 
can be shown using once again the Cardy-Peschel formula(50j, with two angles 2n. |f] This yields another 
logarithmic contribution (c/4) logL. Performing the analytic continuation r — ^ ivpt + e, the result can be 
written as 



-^(0^7 log 



L TTVpt 

— sm — - — 

TT L 



est 



(i»e). 



(33) 



3.3.2. Entanglement entropy We will here compute the EE for the case of an initial cut in the middle of 
the chain La = Lb — L/2. The bipartition coincides with the cut for now. As emphasized before [2], the 
calculation of the EE requires the analytic inversion of the conformal mapping (Eqs. 27 and 28 with r — > 2e 
as in the previous section and the same goes for Fig. [6| . This is important in order for the final analytic 
continuation 9 — > ivpt to be meaningful. In this simple geometry it is indeed possible: 



z{w) = coth 



/ 7re 



1 + C 



\2LJ 1 - C 



with 



c 



/ sinh f (w + e) 
sinh "I (w — e) 



The one point function of the Calabrese-Cardy twist operator must be evaluated at some point w — 
9 real and ~e < 9 < e 



1 dz 



Qm{z) dw 



oc 



TTsinh [f 2e] 



2 sinh [^{e + 9)] sinh [f (e - 6i)] 



(34) 
9, with 

(35) 
(36) 



f One can also use the method exposed here using an infinitesimal transformation L L + SL. 
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Taking the analytic continuation 



— > ivpt, we find that in the limit vpt^e 



L 



+ est 



(37) 



This precise result had already been guessed from numerical simulations in the quantum Ising chain in 
transverse field, see Ref. [5S1I3D]. It is also possible to treat the more general case when the bipartition {i.e. 
the position of the twist operator) does not coincide with the position of the quench, as in [2] for the finite 
case. In our setup the quench still happens at L/2, but the EE is measured for a subsystem A of length 
La = L/2 — I. The twist operator must therefore be evaluated a,t w — 9 + with e < 6 < e. The calculation 
is more cumbersome, but can nevertheless be carried to the end. We find after the analytic continuation 
9 — > it, setting vp = I for convenience. 



S{t) 



12 



log 



^ TT ) 



COS 



2 TTjt+ie) 



COS' 



t) (cos 



2irl 



cosh 



2-JTf 



COS 



27rt 



COS 



2 Tr{t — ie) 



COS 



2 irl 



For t 3> e this expression can be further simplified (we reintroduce vp in the formulae). 

log L + est , < vpt <l or L — l<vpt<L 

S{t) 



\ L 



■)]} 



(38) 



(39) 



est' 



/ < Vpt < L-l 



S{t) is of course periodic with period L/vp. Notice also that the EE is constant at the beginning, until time 
l/vp. A limit of interest is that of an infinite system l,t ^ L: 

{I log I + cst , t < I 

(40) 
f log {t^ - P) + est' , t>l 

where we recover the Calabrese-Cardy[3j result. Actually, one finds here that our results are nothing but 
the Calabrese-Cardy results in an infinite geometry, where vpt and I have been replaced by their "chord 
counterparts" ^ sin (^^-77^) and ^ sin (77). This simple relation with the infinite case is standard in CFT[2 , 
but breaks down as soon as the quench is not located exactly in the middle of the strip (like in Fig. [s] if 
La = £ ^ L/2). Sec Eqs. 51 and 52 for an example where this trick does not apply (with La = L/3, 
Lb = 2i/3). 



S.4- Semi-analytic method for the Loschmidt echo in other geometries 

We are now interested in geometries where the cut is not in the middle of the chain {La ^ Lb). Usually, 
the loss of the left-right (when La = Lb) symmetry makes the problem much more complicated, and a full 
explicit solution for the EE and the LLE becomes less likely. However, we wish to show here that the difficulty 
can be circumvented, at least for the Loschmidt echo. The idea is to perform the analytic continuation semi- 
numerically. We explain this procedure in the case the semi-infinite geometry i — La ^ Lb, L, and will use 
it extensively in the rest of this paper. 



■'I 


1 




T 


• 1 




-6-10 lb 



Figure 7. Semi-infinite geometry (left) and upper lialf-plane ] 
is given by Eg. |41| 



' (right). The conformal transformation w{z) 
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We consider a subsystem of size La = coupled at time t = to a semi-infinite chain. In imaginary 
time, we need the semi-infinite geometry shown in Fig. [Tj The Schwarz-Christoffel transformation mapping 
the upper half-plane to this semi-infinite geometry reads in this case 



'w{z) 



2b^ 1 

1-622 



log 



(f 









(41) 



Our parametrization of w{z) ensures that w'{±b) = 0. The distance r between the two slits is given by 



2e 

T(h) = w{b) — w(—b) = — 

TT 



2b 



1-52 



log 



b-1 



and the variation of the (logarithmic) Loschmidt echo J- 

f{z) = S{w, z){dw / dz)~^ dX z = b: 

(JJ" = ^ Rcs[/(z);z = b] 5t 
= ^Res[/(z);z = 6] ^56 



We get 



F{b) 





2 




I-&2 



log|6|-log|62 



(42) 

log C can be expressed as a residue of 

(43) 
(44) 

(45) 



where 5Re[a] is the real part of a. This is important because we now want to go back in real time. As already 
suggested, one important difference with the two preceding cases is that the equation (42) relating t to 6 



cannot be explicitly inverted. The analytical continuation (r — > ivpt + e) can nevertheless be performed 
numerically. For each time t, we solve 



ivpt 



rib) 



21 

TT 



2b 



I-&2 



log 



1 



6-1 



(46) 



and the solution b(t) can be injected in Eq. 45 Some further technicalities are presented in Appendix C 
3.5. Extension to arbitrary aspect ratios 



3.5.1. Loschmidt echo The method developed in Sec. 3.4 can be pushed further, to treat any aspect ratio 
X = (./ L = La/ {La + Lb). The conformal transformation which maps the upper half plane onto the 
"asymmetric double-pants" geometry in Fig.|8]is given by 



w{z) 



L 

TT 



a; log 



z - 1 



+ (1 - x) log 



z 



a 



z + a 



(47) 





jr ( 


L-i 






—a —b —I 



1 b a 



Figure 8. "Asymmetric double-pants" geometry (left) and upper half-plane H (right). The conformal 
transformation w{z) is given by Eq. |47| 
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As before, we first need to determine the distance between the two sUts t = w{b) — w{—b), where b is the 
solution of w'{ztb) — 0. We find 



(1 — x)a + X 



2L 



a; log 



b-1 
bTl 



+ (1 - a;) log 



,-b 



(48) 
(49) 



Then the free energy can again be computed as a residue at z 



T{a) = 



6 



log \a\ 



3x^ - 3a; + 1 
a:(l — x) 



log|a + l|- 



b, and we get after some algebra 
a;^ — a; + 1 



a;(l — x) 



log|a-l| 



+ ^ log \ax + (1 - a;)| + ^ log |a(l - a;) + a:| 



(50) 



When X is the inverse of some integer {x = 1/p ,p G N), getting a as a function of exp {7tt/2L) amounts to 
finding the roots of a polynomial of degree 2p — 2. Doing so for r = ivpt + e we get a whole set of solutions, 
but only one of them corresponds to the actual real-time LLE. 



3.5.2. Entanglement entropy The entanglement entropy turns out to be much more complicated in the 
asymmetric geometries. As already stated, the calculation of the one-point function requires the inversion of 



the conformal mapping of Eq. 47 When x = 1/q where q is some integer, inverting the mapping amounts to 
finding the roots of a polynomial of degree q. An analytic expression of those is in principle available only for 
q = 3 and (7 = 4. Such an analytic expression is needed in order to perform the analytic continuation to real 
time. In practice, the expressions we get for the entanglement entropy are gruesome, and not particularly 
illuminating. We managed to treat the a; = 1/3 case, and present our result here for the sake of completeness. 
First we introduce the auxiliary function 

|4(1 - 3e) [17 - cos(37rt)] cos (2f^) - 24^^2(1 - e) ^1 + 16e - cos(37rt) sin (^f^) |^^^ 

(4 + 2e) sm (^) 

We denote by = 5ie(/e) and — 5m(/c) the real and imaginary parts of /g. We also set vpL = 1 for 
convenience. The EE can then be expressed using f^{t) and its derivative f^{t) — g'^{t) + ih[{t): 



c 



[Kit)]' 



S{t) = —log { ' " ^ ^Kcst- (e<t<2/3), (52) 



2 



where cscu — 1/sinu. It is important to keep the regulator e finite in Eq. 52 because the denominator is of 
order e. We also note that this expression can probably be further simplified. A comparison with numerical 
simulations is shown in Sec. ID 



4. Light-cone effects: plateaus of the LLE and non-universal decay of the entanglement 
entropy 

We wish now to interpret some of the results we derived for the Loschmidt echo, as well as for the 
entanglement entropy. As emphasized in the introduction, we are focusing here mainly on finite-size effects. 
This is especially interesting in the case of a local quench, because the quasiparticles emitted at t = can 
bounce back on the boundaries. Now something interesting happens when the left-right periodicity is broken 
La 7^ Lb: for certain periods of time, all the quasiparticles can be in the same subsystem (see Fig.jo]). Since 
they are supposed to transport information about entanglement and correlations, it is reasonable to expect 
that the behaviour of both the (logarithmic) Loschmidt echo and the EE will drastically change. This is 
indeed the case, both quantities exhibit light-cone effects. However, we will see that they are quite different. 
For the EE we see a slow non symmetric decay, whereas we get striking plateaus for the LLE. 
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Figure 9. Drawing of quasi-particles emitted at t = by the local quench. They propagate at the speed of 
light vp and bounce against the wall. Left : x = 1/2, centre : x = 1/3, right : x = 1/5. The plateau regions 
(filled in grey) correspond to the case when the two quasiparticles are in the same subsystem. 



We also wish to compare our analytic results to several numerical simulations in free fermionic spin chains, 
which are described in the continuum limit by CFT with central charges c=lorc=l/2. This section is 
organized as follows. First we give some details about the model - so called XY-chain in transverse field - and 
show how the quasiparticles can be identified. We then discuss separately the behaviour of the entanglement 
entropy and the Loschmidt echo. Finally we introduce a variant of the Loschmidt echo ( "detector geometry" ) 
specifically constructed to probe the light-cone effect, and discuss the results. 



XY spin chain in transverse field 
We consider the following Hamiltonian 

^-,t-h (53) 
i=i 

where the boundary conditions for the spins are free and L is assumed to be even. As is well-known, H 
can also be expressed as a quadratic form in fermionic operators, which makes it one of the simplest exactly 
solvable spin chains. The fermionic structure also allows to calculate numerically both the entanglement 
entropy and the Loschmidt echo for very large system sizes. The phase diagram is quite rich. In particular 
there are two critical lines, lying in different universality classes. 



H 



11 y 



J^.l.l. XX chain: free boson critical line This corresponds to r = and — 1 < /i < 1. H can in this case 
be rewritten up to a unimportant constant as 

-ff = ^ e™4„c™ , e„i = -h- cos f j , (54) 

m=l \ + / 

where we have the usual anticommutation relations {cm,cjj} — Smn and {cm,c„} = {c}„,c}j} = for the 
fermions. The quasi-momenta are the fc,„ — mir /{L -\- 1) and the Brillouin zone extends over the interval 
(0, tt). At zero magnetic field /i = 0, the ground-state can simply be obtained by filling half of the first 
quasi-momenta: 

E,{L) = - E cos (^) (55) 

The asymptotics of when L — > oo can be accessed using the Euler-Maclaurin formula: 

E^{L)^aL + P-^+0{llL^) (56) 

a and /3 are non universal constants and direct comparison with Eq. [3] yields vp y-c= 1. The Fermi velocity 
can simply be calculated from the dispersion relation, vp = \ ^\\k=k =,r/2 ~ ^ ^'^'^ S®* textbook result 
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that the XX chain lies in the c — 1 (free boson) universahty class. As already explained, we expect that the 
local quench will mainly generate low-energy excitations, therefore at energies close to the Fermi level. In 
this region, the dispersion relation can be linearised, and this is precisely what the CFT describes. 



At non-zero magnetic field, the filling fraction is not 1 /2 anymore and the number of fermions is an integer 
n 7^ L/2. This case turns out to be trickier. Let us focus for the moment on the case when n = pL, where p 
is some rational number. The ground-state energy is 

Eo{L) = - £ cos (^) , (57) 

m=l \ + / 

and applying the Euler-Maclaurin formula we get (vp = sin kp ^ sinpTr) 

Eo{L) = a'L + p'~^[l- I2{p l/2f] + 0{l/L^) (58) 

To understand this result and its implications, we follow Zagoskin and Affleck [5T]. Let us first define 
our convention for the Fermi level. We choose to place it between the last occupied and the first empty 
wave- vector 

n + 1/2 pL + 1/2 

kp — TT = TT (59) 

L+l L+l ^ ' 

Except at half-filling p = 1/2, the particle-hole symmetry is not exact on the lattice: e(fc„) — e{kp) ^ 
£{kp) — e{kn+i). This effect has surprisingly important consequences in the continuum limit. In a bosonisation 
picture, it acts as a phase shift |51| and changes the effective boundary conditions on the boson field. 

As explained in Ref. 51j, this difficulty can usually be circumvented by choosing appropriate conventions 
for the length of the chain and the number n of fermions. For example at quarter-filling, choosing the length 
of the chain as L = 4n -|- 1 restores the particle-hole symmetry and therefore cancels the anomalous term 
12 X (1/4 — 1/2)^. This yields the desirable result c = 1 for the central charge. In this case the boundary 
conditions at both ends of the chain are the same, which is easier to treat from a CFT perspective. However, 
we are interested here in the effect of local quench, where the interaction between two subsystems A and B 
is switched on. For parity reasons this trick cannot be applied to the chains A, B and A U i? at the same 
time. In other words, it is not possible (away from half-filling) to maintain the same boundary condition 
everywhere and we must introduce boundary changing operators in the CFT calculation. It is known that 
they don't affect the leading term in the entanglement entropy, and we therefore expect our formulas to 
hold away from half-filling in that case. However, these operators can have a noticeable effect on the scaling 
behaviour of the (logarithmic) Loschmidt echo and we expect our results to break down away from half- 
filling. Although these effects can be treated by standard CFT techniques, we do not study them in this 
paper. Therefore, we will stick to numerical checks at half-filling. Finally, let us mention that this discussion 
also applies to the XXZ chain in a magnetic field. There is also a phase shift away from half-filling. It takes 
a more complicated form which has been computed using Bethe-ansatz techniques |52l [55] . 

4-.1.2. Ising critical line This corresponds to /i = 1 and r > 0. In this case the Hamiltonian can be 
diagonalized by a Bogoliubov transformation: 

L 

- ^ {dld^ - 1/2) (60) 

m— 1 

and the one-particle energies are given by 

Em = \/ih + COS kmf + sin^ fc„ (61) 

The wave- vectors km can be obtained as solutions of an implicit equation, and lie in the interval (0, tt). For 
r = \ we have the Ising chain in transverse field, and they take a simple form: 

771 TT 
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This model is known to have central charge c = 1/2 for all values of r > 0. In the continuum limit, the 
dispersion relation is 



e(fc) = y^(l + COS fc)2 + si 



sin^ fc, 



(63) 



II 



the Fermi momentum is at kp = tt, the Fermi velocity is vp = de/dk\k=TT = r, and the CFT again describes 
excitations close to the Fermi point. For r < a/3/2, something interesting happens: e{k) is not convex and 
the Fermi speed is not the fastest speed available for the excitations (see Fig. [lOjfor an example). This 
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Figure 10. Left panel: group velocity v{k) = \de/dk\ for three different values of r (r = 0.6,1,1.4). 
fmax is the maximum group velocity allowed. Right panel: ratio Vmax/vp, where vp is the Fermi velocity 
Vp = v{k = n). For r < \/3/2, the ratio is greater than one. 



maximum speed 



~ max v(k) = max 

k k 



de 
dk 



(64) 



is sometimes called the Lieb-Robinson speed [Hj and corresponds to the maximum velocity allowed for 
the propagation of a signal in the quantum system. Its role in the context of global quenches has been 
emphasized[55l[56l[57]. As we shall see, these fast particles — which propagate faster than the CFT particles 
discussed in the introduction — affect both the entanglement entropy and the Loschmidt echo. Conversely, 
in the case r > V3/2, the Fermi speed is the fastest available speed, and the behaviour of both the EE and 
the Loschmidt echo is rather affected by particles which propagate slowly (compared to vp). 



4-2. Non-universal decay of the entanglement entropy 

We first focus our attention on the EE. The only asymmetric finite-size system we have been able to 
treat corresponds to the aspect ratio x = La/L = i/L = 1/3. The solution given by Eqs. 51 and 52 can 
be compared to numerical simulations in free fermionic systems. Fig. [TT] shows such a comparison in the 
XX chain for very large system sizes, up to L = 3072. The agreement between the two is excellent at the 
beginning, but suddenly breaks down precisely at the time {t = 2La/vp) when the two light-cone boundaries 
(or position of left-moving and right-moving CFT quasi-particles) lie both in subsystem B (see Fig.|9]in the 
middle). Then the entropy slowly decays up to time t = 2Lb/vp, when the CFT calculation becomes 
accurate again. Such a decay was already observed in the semi-infinite geometry in Ref. [51], also for the 
XX chain. It can be attributed to the slower quasi-particles progressively entering the subsystem B. 



Several interesting observations can be made: 

• Global scaling with the system size L: in the "CFT region" (corresponding to t < 2La/vp and 
t > 2Lb/vp in the first period) the numerical simulations show a global scaling (c/3)logL. This is 
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Figure 11. Entanglement entropy for the aspect ratio x = 1/3. Left panel : Numerical simulations 
in the XX chain, shown against the rescaled time {vp/L)t. Right panel : zoom on the "CFT region" 
< t < 2L/{3vp), and comparison with the analytical calculation of Eqs. |51|and|52| 



a consequence of the presence of two slits in the CFT setup, which contribute bothjl] to a (c/6) log L for 
open chains. This is not true in the "decay region" {2La/vp < t < 2Lb/vp) where the global scaling 
is half of the previous one: (c/6) log L. This can be interpreted as follows. From the point of view of 
the subsystem A, there are no CFT quasi-particles anymore in A and thus A sees — up to non-universal 
contributions — a wave- function close to the ground-state wave function \AU B) of the total system, 
which contributes exactly to a (c/6) log L. 

At t = 2Lb, the CFT particles are again in their own subsystem and contribute to the EE as predicted 
by CFT. However, one could expect non universal contributions to deteriorate the agreement with the 
numerics. This effect seems to be negligible (Fig. [TTj). 



To test the universality of the EE, and the effect of the non-CFT quasiparticles, we also performed numerical 
checks in the XY chain in transverse field, for various values of the r parameter in the Ising critical line. As 
explained in the previous section, the case r < a/3/2 is particularly interesting, and puts the CFT calculation 



to a more serious test. Some results are presented in Fig. 12 for r = 0.6 and the aspects ratios x = 1/2 and 
X = 1/3. 



As explained in Sec. |4.1.2[ the maximum speed is not the Fermi speed in this case, and some quasi 
particles will therefore enter the subsystem B before their CFT counterpart do. Although the effect of these 
fast particles should be small, it does not seem to disappear in the thermodynamic limit. For example in 
the symmetric case x = 1/2, there is a priori no decay region, and one expects the EE to be periodic with 
period T — L/vp. However, at time t — i/wmax the fast particles suddenly change subsystem, and a small 
deflagration — analogous to the one observed in [57] for a global quench — can clearly be seen on the EE. The 
CFT calculation does not seem to hold anymore, due to these non universal effects. The same phenomenon 
happens at time t — 2L/vma^- For the aspect ratio x = 1/3 the same interpretation holds. Their influence 
is however quite visible in the decay region. 



4-. 3. Plateaus oj the LLE 



We now turn to the (logarithmic) Loschmidt echo. In the symmetric case La — Lb, its behaviour is very 
similar to the EE, and the two quantities are given by a very similar formula (33)-(37|. For an asymmetric 
system La ^ Lb such a close similarity between the EE and the LLE does not hold anymore, especially when 
the two quasiparticles are in the same subsystem {i.e. the "decay region"). The behaviour of the LLE is 
particularly striking: we get very flat plateaus, which can be accurately calculated in CFT (Eqs. [48| and [50| . 
See for example Fig. 13 for a comparison with the numerics in the cases x = £/L = La/{La + Lb) = 1/3 



17 




Figure 12. Left panel: Entanglement entropy as a function of the rescaled time T = {vp/L)t for the 
ratio X = 1/2, the anisotropy r = 0.6 and for systems sizes L = 128,256,512,1024. At rescaled times 
fa = vp/vma,^ and 2Ta the CFT calculation breaks down, although it remains a good approximation. Right 
panel: Entanglement entropy for x = 1/3, r = 0.6 and for system sizes L = 192, 384, 768, 1536. The rescaled 
times of interest are Ti = 2t)i?/(3i)max), 2ti and 3ti. 
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Figure 13. Logarithmic Loschmidt echo J^{t) = —\ogC{t) for the aspect ratios x = La/{La + Lg) = 1/3 
(left panel) and a; = 1/5 (right panel). Red curve : CFT calculation. Black points : numerical calculation 
in the XX chain. 



and X = 1/5. Contrary to the EE, these plateaus are also symmetric with respect to t i — > 2L/vp — t, and 
have the same global scaling (c/4)logL as the other parts of the curve. To check the universality of our 
calculation, we also performed numerical checks in the XY chain in the case Wmax > vp- The results are 
shown in Figure 14 for the aspect ratios x = 1/2 and x 



: for the aspect ratios x = 1/2 and x = 1/3. Interestingly, the Lieb- Robinson quasiparticles 
have a different effect on the LLE, compared to the EE. Finite size effects are enhanced when they change 
subsystem, but these diminish when increasing L. For very big system sizes the numerical calculations match 
very well our analytical prediction, even in the "plateau region" . It therefore seems that the Loschmidt echo 
is less sensitive to non-universal effects than the EE. 



It is also interesting to study the semi-infinite geometry La ^ £ <^ L = La + Lb ■ In this case the period 
is infinite, and one gets a single plateau starting at time t = 2£/vp. Although the analytical solution is given 
in implicit form, we can extract some interesting quantitative features about the small time and long time 
behaviour of the LLE. J^{t) — — log£(t) is shown in Fig. 15 We discuss now this curve (Eqs. 45 and 461 in 
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{vFlL)t = x {VF/L)t = x 



Figure 14. Left panel: Logarithmic Loschmidt echo as a function of rescaled time T = (vp/L)t for the 
aspect ratio x = 1/2, the anisotropy r = 0.6 and for system sizes L = 256,512,1024,2048. Right paneh 
Logarithmic Loschmidt echo as a function of rescaled time T = {vp/L)t for the aspect ratio x = 1/3, the 
anisotropy r = 0.6 and for system sizes L = 384, 768, 1536, 3072. In both cases, finite size effects decrease 
for bigger and bigger system sizes. Note that, contrary to the entanglement entropy, the presence of fast 
particles (i)max > "f) does not spoil the plateaus. 




{vFlL)t = X {VF/L)t = x 



Figure 15. Left panel: (logarithmic) Loschmidt echo T{t) = — \ogC{t) for the semi-infinite geometry and 
central charge c = 1. Right panel: zoom on the plateau, for vpt/L > 2. The red curve is Tit) and the thin 
blue curve is the asymptotic expansion of Eq. |67| 



some detail. At small time t <^ £, we recover the expected behaviour [12] 

J'(t) = ^logi + cst. (65) 

This is simply a consequence of the fact that the quasiparticles do not feel the left boundary yet, and behave 
as if the system was infinite. Then J^{t) keeps growing, up to a certain time when it reaches a maximum. 
This time tmax can even be calculated exactly 

— ~ 1.3538 X — (66) 
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Also, the long time behaviour can be extracted [§] 

^^'^ - ^^^^ ^ 6nHvlt/er + 

The decay is algebraic, with a small — but universal — prefactor. This is actually the main reason why the 
plateau is so flat. 

4-.4- "Detector" geometries 

To further probe the light-cone effects, and confirm the appearance of the plateaus as a generic phenomenon, 
we introduce here a variant of the (logarithmic) Loschmidt echo. We call it "detector" for reasons which will 
become transparent later. It is defined as 

V{t) = - log\{A(E) B\e^"^\C (gi D (I) E)f (68) 
The difference is that the wave function after a time t I'lpit)) = e*^*|A(g) B) is not compared with the initial 



state, but with another tensor product of 3 ground-states |C), \D), and \E), see Fig. 16 For simplicity we 
only consider the case La = Lb = L/2 and Lc = Le here, so that the only remaining free length in the 
problem is the middle one Lo ■ Now what should be the behaviour of such a quantity? Obviously the quench 

Final state ^^^~^^^^^^^^~^^^^^^^^~^^^H 



Lc — Le Lo Le — Lc 



Initial state ^^^^^^"^^^^^^^♦^^^^^^"^^^^^^H 



La = L/2 Lb = L/2 



Figure 16. The system is initially prepared in the ground state \A (g) B) of Hy^ + Hg, then at time t > 
it evolves with the total Hamiltonian H = Hj^ + Hg -\- H^g. Finally it is compared with the ground state 
|C ^5 (gi i?) of He + H]j + He- Blue dashed arrows are a graphical representation of the quasiparticles 
emitted after the quench. 

at < = is the same, so that there are again quasi-particles emitted at the interface between subsystem A 
and B. They propagate at speed vp, but from the point of view of the final state \C ® D ® E)^ they are at 
the beginning in the same subsystem, namely D. This situation is particularly interesting, because it means 
that the detector 'D{t) should start with a plateau, up to a time t — i£)/2, when the quasi particles emitted 
to the left and to the right simultaneously get out of D. In this sense, the final state is able to "detect" the 
quasiparticles. One other advantage is that I'(i) can also be computed in CFT, using the same techniques 
as for the Loschmidt echo. The technical details are presented in [Appendix B[ We present here only the 
results, focusing on two geometries: 

• The "symmetric detector" Lo = L/2 (and therefore Lc — Le ^ L/A). In this case the enhanced 
symmetry even allows a full explicit solution (we set vp = 1): 



V{t) - - 



ilog|2cos27rt/i| -f 21og 1 - e'"'^' ^ ^j2cos2T:t/ L 



,0<t<L/2 (69) 



This happens to be the only geometry in which the plateau can be obtained in explicit form. V also 
satisfies 'D{L — t) = I?(i), I?(t + L) = Vit) and is quadratic at small times I?(t) oc i^, similar to the 
Loschmidt echo at long times for the semi-infinite geometry. 

§ For t — ^ oo, the solution a{t) of Eg. |46| can be written as a{t) = ^^lyrj, where b{t) = a(t) — i/3{t) is a complex number with 
very small modulus. One can then obtain the following asymptotic expansions: /3{t) ~ — t — ^, . . and a(t) ~ — 2/3(t)^ log/9(t). 
Injecting this in Eq. 1451 and expanding around a = 0, we get the result. 
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• The "infinite detector" ^ £ L. The solution is this time not fully explicit, but 'D{t) can 
nevertheless be expressed using the Lambert W function, solution in the complex plane of the implicit 
equation tee™ — z : 

V{t) ^ Uog\a{t)[ait) ~ l]\ , a(<) - ly (e-2*/^-i) + 1 (70) 

As for the Loschmidt echo in the semi-infinite geometry, several interesting features can be extracted 
from this implicit solution. At small times 'D{t) behaves as oc i^, and the plateau lasts up to i = £/2, 
consistent with the quasiparticlc interpretation. After this, the quasiparticles (and 2?(<)) do not feel 
anymore the subtleties of the detector geometry, and we gel[j||for t/£ ^ 1: 

V{t)^Uogt. (71) 

We recover the result of the infinite geometry. 

We also performed numerical tests of the formula Eqs . [69| and [70| in the XX chain and for very large system 
sizes. The numerical results in Fig. |17| are in excellent agreement with our analytical calculations. 




{vF/L)t = T (vF/e)t = T 

Figure 17. Left panel: symmetric detector. Black dots are the numerical points for system sizes 
L = 1024, 2048, 4096, whereas red curves are the CFT prediction of Eq. |69| Right panel: infinite detector 
for I = 256 and L = Lao = 4096 S> Black dots are the numerical points, red curve is the prediction of 
Eq.[70] 



5. Conclusion 

In this paper, we study the time evolution of a one-dimensional quantum critical system after a local quench. 
This "cut and glue^^ quench is particular, in the sense that it corresponds to a single conformal boundary 
condition in the scaling limit. This is clearly restrictive, and the problem of quenches corresponding to 
different conformal boundary conditions (and more generically to the possible RG flows between those) 
deserves some further investigation. For such a quench, we have argued that the evolution of the system is 
periodic, with period 2L/vf if L is the size of the system and vp is the Fermi velocity. This is true only in 
the continuum limit, but we observe that it is a very good approximation for large finite-size systems, which 
typically exhibit a number oscillations just after the quench. For very long times however (t ^ L/vp), the 
finite-size effects spoil the periodicity of the system and the CFT description is not accurate anymore. 

The restriction to this particular type of quench also allows us to get simple universal results which 
complete the ones obtained over the last years by Calabrese and Cardy [H |2] and Eisler and Peschel [3] . 
Moreover, we aimed to give some further evidence that quantities defined as simple overlaps can mimic 
the behaviour of the EE [TH]. Our motivation comes partly from the fact that such overlaps should be, in 

II The solution a(t) = —a{t) + i(S{t) admits the foUowring asymptotic expansion for t ^ £: a{t) ~ log{TTt/£) and /3(t) ~ TTt/£. 
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principle, much easier to measure in real experiments than the EE or Renyi entropies. Indeed, such overlaps 
are nothing but the probability of observing the system in its ground state after a quench H ^ H' . Although 
Cardy [18] recently proposed a type of local quench that would lead to the measure of the Renyi entropies, 
it seems to us that the "cut and glue" quench is still more likely to lead to a realistic implementation. In 
this paper we thus focus on overlaps, and more precisely on the (logarithmic) Loschmidt echo (which is a 
time-dependent version of the LBF introduced in |19|). and compare its behaviour with the one of the EE. 
From a CFT perspective, both quantities can be computed through analytic continuation r — vpt+ie, where 
e is an UV cutoff. However, for general geometries, it is not possible to perform this analytic continuation 
to compute the EE (because of the absence of an analytic expression for the one-point function in generic 
cases). On the contrary, we find that for the LLE, a semi-analytic treatment involving a numerical solution 
gives very accurate results. 

We provide numerical checks of our analytic calculations in free fermion systems. We observe light- 
cone effects, both for the EE and for the LLE. These effects, however, are not quite the same for the two 
quantities. For a "cut and glue" quench Ha + Hb H = Ha + Hb H- ^^ab {^a < Lb), within the first 
oscillation of the system < t < 2L/vp, one can distinguish two regions. The "CFT region" corresponds to 
t < 2La/vf or t > 2Lb/vf, where we observe that both the EE and the LLE are very well approximated 
by the CFT calculation. Then the "decay region" or "plateau region" {2La/vf < t < 2Lb/vf), where the 
discrepancy between the behaviour of the EE and the LLE is more important. The EE decays slowly, as 
observed first by Eisler and Peschel j3j in a semi-infinite system, who interpreted this decay as an effect 
due to slow particles. The LLE does not exhibit such a decay, but rather stays almost constant, leading 
to pictures with flat plateaus. For large systems these plateaus are very well approximated by the CFT 
result, which predicts a behaviour that is symmetric under t i— ^ 2L/vf — t. The EE is more sensitive to 
non-universal effects than the LLE for this "cut and glue" quench. We observe, however, that the Loschmidt 
echo depends on the filling fraction of the fermion systems which we investigate in the c = 1 case. The CFT 
results we have obtained in this paper for the LLE turn out to be correct only at half-filling, and these must 



be modified for other filling fractions, as discussed in section 4.1.1 (the same phenomenon appeared already 
for the LBF). 

Finally, it would be interesting to have checks of all the CFT results — which have accumulated over the 
past years — coming from other fields and using different techniques such as integrable systems [53 [59l [60] 
or numerical calculations in interacting systems. Another important direction is to study how these CFT 
results can be related to more general types of quenches, where the parameters in the Hamiltonian are not 
turned on/off instantaneously, but rather continuously on a short time scale. It is known that the behaviour 
of the system can be affected in different ways at short time (see for example [13]). CFT techniques lead 
to results that are valid only on longer time scales. How to relate these different regimes is an interesting 
open question. More generally, it is not clear to us what are the most relevant physical aspects which are 
missed in these CFT models, and how one could make the latter more realistic. This seems to be of crucial 
importance if one believes that such results are to be compared with the ones of an actual experiment. 
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Appendix A. Loschmidt echo and free fermions 

Appendix A.l. Diagonalization 

The general method to compute of the EE is known since Peschel [61j . For the type of local quenches we 
are interested in, we refer to Refs.[3|, |4T]. where the procedure is explained in details. Calculating the EE 
usually amounts to diagonalizing a 2L^ x 2La matrix. 
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We give here some details about the numerical computation of the Loschmidt echo in free fermionic 
systems. We consider the XY chain in transverse field 

A Jordan- Wigner transformation 

a] + ig] = 2c] exp c^q^ (A.2) 

a] = 2c]cj - 1 (A.3) 

then allows to express the Hamiltonian as a quadratic form in the fermionic operators 

which can be diagonalized using a Bogoliubov transformation: 

L ^ 
/m — E/ '^O'^'^j ~^ "^jmCj ' /m ~ (/m) (^-5) 

i=i 

The Hamiltonian is then rewritten as 

HauB =Ea+J2 ^mflfra. (A.6) 
m 

where we can choose all the one-particle eigenenergies to be positive > Vm. Eq is therefore the ground- 
state energy. Denoting by |0) the fermion vacuum (cj|0) = 0), the ground-state of the total chain can then 
be expressed as 

\AyjB)= hh...fL\0). (A.7) 

One can also do the same for Ha^b- Denoting by {c?m,c?J,i} the set of fermions which diagonalize the 
Hamiltonian, 

the ground-state is 

\A(g)B) =did2...dL\0). (A.9) 

Appendix A.2. Loschmidt echo and bipartite fidelity 
Recall that the Loschmidt echo is given by 

Cit) = \{A(x)B\exp{itHAuB)\A(x)B)f (A.IO) 
To compute it, it suffice to notice that the state 

\n) =di...dL\A[JB) (A.ll) 

is an eigenstate of Ha^b with the minimum eigenenergy Sq. Since the ground state of such an Hamiltonian 
is unique (Perron Frobenius theorem), is actually proportional to \ A® B). We therefore have 

{A U B\d\ . . . d\e^*"^^^di . . . dL\A U B) ^ 

C{t) = ^ (A.12) 

{A\JB\d\...d\di...dL\A\JB) 

The denominator is nothing but \{A U B\A ® the (squared) bipartite fidelity[T9], and ensures 

£(^ = 0) — 1. The numerator N{t) can be expressed using the fermions operators 

N{t) = \ {d\...d\d^{t)...dL{t))\^ (A.13) 
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where d,„(t) = e^'-^^^^rfme'^'^-^us a^d we have used the fact that e'*^-^uB|^y 5^ contributes just to a 
phase. This correlator may be reduced to a Pfafhan using Wick's theorem. Using the relation Pf^ — det, 



det 



(A.14) 



where M is antisymmetric, while iV^*^ is symmetric. The matrix elements are given by 

M,i = {d]d\) = -{djdi) = -{d,{t)di{t)) (A.15) 

Nf^ = {d]d,{t)) (A.16) 

In the determinant, adding the second column to the first and then the first row to the second we find 

7V(*) + M iV(*) 
iV(*) - M 



U{t) = 



dot 



det(iV(*) + M) 



(A.17) 



Appendix A. 3. Bipartite fidelity 



The bipartite fidelitvP^ is the (squared) overlap \{A\J B\A® B)]"^ , and is a direct by-product of the preceding 
calculation. It can therefore be expressed as L x L determinant: 



\{AuB\A(g)B) 



det 



(A.18) 



This result was actually already used in Ref. p[H]- 



Appendix B. Details on the detector geometries 

Appendix B.l. The "symmetric detector" 

This turns out to be simplest detector geometry. The conformal transformation which maps the upper 



half-plane onto a strip with 3 slits of Fig. Bl 



w{z) = ^ [log(l - z^) - 21og(l - a^z^)] (B.l) 




Figure Bl. "Detector" geometry (left) and upper half plane H (right). The conformal mapping from the 
latter to the former is given by Eq. |B.1| 



This mapping has singularities ai z — ±1, z = ±a and extrema at z = 0, z = ±-\/2 — l/a^. The distance r 
between the two slits is 

r(a) = Aiog(4a2(a2-l)) (B.2) 

The crucial point is that this bisquare relation can easily be inverted. We then use an infinitesimal 
transformation w 1-^ {1 + 5t/t)w at the right of path T, w 1— > w at the left. Denoting as before by 
/(z) — S{w, z){d'w/dz)^^ , the variation of the free energy is 

^ = f / fi^)dz (B.3) 
or DTT Jr' 
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where T' is the image of T under the inverse map w i-^ z (see Fig. Bl). The curve T' can be closed onto a 
contour using its mirror image along the real axis, and the calculation boils down to a residue. 

SV 



- - xRcs[/(z);z = 0], 

OT D 



SO that 



loe 



(a2 - l)(2a2 _ 1) 



Inverting Eq. |B.2| with t ~ it wc finally get 



-log|2cos27rt/i| +21og 1 - e"*/-^ V2cos27ri/^ 



,0<t<L/2 



(B.4) 
(B.5) 

(B.6) 



and D{L — t) = 'D{t + L) = !?(<). In this equation the square root has to be understood as y/a = i^J—a for 
real negative a. 



Appendix B.2. The "infinite detector" 

Here the required conformal transformation reads 

i r a- 1 



w{z) 



2-K 1 



- log (z - 1) - log(z + 1) 



(B.7) 




Figure B2. "Infinite Detector" geometry (left) and upper half plane H (right). The conformal mapping 
from the latter to the former is given by Eq. |B.7| 



has singularities at z — ±1 and extrema at z = 0, z = ±y^. The distance between the two slits is 

r(a) = ^[a + \og{a - 1)] (B.8) 
The detector in imaginary time is solution of the equation 

— - ^ xRes[/(z);z = 0], (B.9) 
which can easily be integrated 

I?(a) = ^log|a(a-l)| (B.IO) 



Contrary to the previous case, Eq. |B.8| cannot be explicitly inverted. The solution can nevertheless be 
expressed as 

a(t) = 1 + W^(e"2*/^-i) (B.ll) 

where W is the (multivalued) Lambert W function, solution of the implicit equation we^ — z. This is 
especially useful when trying to evaluate the Detector, because W is part of the toolbox of usual symbolic 
computer programs. 
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Appendix C. Intriguing trajectories in the complex plane 



In all cases we encountered in this paper, the Loschmidt echo (or the detector) could be expressed as 

- log C{t) = - log £ [ait)], (C.l) 

where except in favourable cases a{t) is the solution of some implicit equation in the complex plane. The 
goal of this section is to make this solution more intuitive, by simply showing this parametric function in 
the complex plane, especially at times corresponding to a plateau. 



We focus on three geometries : the asymmetric Loschmidt echo for x = £/L = 1/3, the "symmetric 



detector", and finally the "infinite detector". The results are summarized in Fig. CI where a{t) is shown 



for these three geometries. For the Loschmidt echo we also show the curve solution of C{a) — C [a{t — L)], 
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Figure CI. Left panel: Loschmidt echo for x = £/L = 1/3. Middle panel: "symmetric detector". Right 
panel: "infinite detector". Black curve is a{t). The fat red curve is Red curve is the iso Loschmidt echo 
C{a) = C [a{t = L)] for the left panel, and the iso Detector ©(a) = T> [a(t = 0)] for the two other. The 
plateau correspond to times when the black curve a(t) almost follows the red one. 



namely the set of a which would give the same values as the observed one in the middle of the plateau 
e.t t — L. This is the red curve in the left panel of Fig. |C1[ Since the plateau region is close to perfect, 
one expects the solution a{t) to be very close to this curve at times 2L^ < t < 2Lb, and this is indeed 
what happens. We also performed the same analysis for the detector, now showing the curves solution 
of I?(a) = I? [(a(t = 0))], since the plateau begins at time t = in this case. We also observe the same 
phenomenon. As an aside, we note that the red "iso-detector" curve is a well-known Cassini oval. 
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